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A particle driven by deterministic chaos and moving in a spatially extended environment can exhibit normal diffusion, 
with its mean square displacement growing proportional to the time. Here we consider the dependence of the diffusion 
coefficient on the size and the position of areas of phase space linking spatial regions ('holes') in a class of simple one- 
dimensional, periodically lifted maps. The parameter dependent diffusion coefficient can be obtained analytically via 
a Taylor-Green-Kubo formula in terms of a functional recursion relation. We find that the diffusion coefficient varies 
non-monotonically with the size of a hole and its position, which implies that a diffusion coefficient can increase by 
making the hole smaller We derive analytic formulas for small holes in terms of periodic orbits covered by the holes. 
The asymptotic regimes that we observe show deviations from the standard stochastic random walk approximation. The 
escape rate of the corresponding open system is also calculated. The resulting parameter dependencies are compared 
with the ones for the diffusion coefficient and explained in terms of periodic orbits. 
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The diffusion coefficient of a pliysical system involves a 
macroscopic measurement of tlie linear increase in the 
mean square displacement of an ensemble of microscopic 
points. The source of this increase is often modeled as the 
result of many random colUsions between the elements of a 
system, for example caused by thermal motions. However, 
on a molecular level, the rules governing the colUsions are 
deterministic. In order to obtain a full understanding of 
the process of diffusion, one must take these determinis- 
tic rules into account. The discovery of chaotic dynamical 
systems which exhibit diffusion has helped facilitate this, 
as one can consider a relatively simple, low-dimensional 
setting in which diffusion can be studied along with the de- 
terministic rules at the heart of diffusion. In this setting, 
one can consider how the diffusion coefficient reacts to the 
variation of control parameters, and hence gain under- 
standing of more compUcated diffusion processes. Here 
we consider a simple, one-dimensional dynamical system 
which consists of an interval map, periodically copied over 
the real line with each interval connected to its neighbors 
via two regions which we call holes. We analytically derive 
the diffusion coefficient for this system and study the ef- 
fects of the size and position of the holes. We show that the 
position of a hole has a crucial impact on the diffusion co- 
efficient, sometimes dominating that of the size of a hole. 
We further study the diffusion coefficient in the asymp- 
totic regime of small holes. We find analytic expressions 
for the diffusion coefficient in terms of the limiting behav- 
ior of a hole, in particular whether they contain a short 
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periodic orbit. This allows construction of a periodic or- 
bit expansion for small but finite holes generalizing simple 
random walk theory, which is based on an effectively un- 
corrected dynamics. Finally, we compare our diffusion 
results with those obtained previously for the escape rate 
and discuss the similarities and differences that arise. 



I. INTRODUCTION 

Recently there has been a surge of interest from mathemati- 
cians and physicists on dynamical systems with holes, that 
is, subsets of phase space that allow trajectories to leak out. 
Varying the position of such holes in simple dynamical sys- 
tems has exposed a strong link between the average lifetime of 
chaotically transient orbits and the location of short periodic 
orbits'. Work on the escape properties of the doubling map 
and related systems has revealed fundamental results for fi- 
nite times, non-monotonic dependence of escape rates on hole 
sizes, and a precise slowing of escape for small holes contain- 
ing a short periodic orbit". The small hole effects have been 
generalized^^, and noise effects have also been considered^-^; 
for reviews see Refs.8-10. This work begs the question: How 
do transport processes such as diffusion depend upon the posi- 
tion and size of a hole? We answer this question for the diffu- 
sion coefficient of an extended dynamical system, interpreting 
the holes now as the links between spatially separate regions 
of phase space. Diffusion is a fundamental transport process 
of many-particle systems, the study of which cross-links trans- 
port theory in statistical mechanics with dynamical systems 
theory'^"^"*. The linear increase in time of the mean square 
displacement of an ensemble of points in a chaotic, determin- 
istic dynamical system is a process known as 'chaotic diffu- 
sion' or 'deterministic diffusion' Chaotic diffusion pro- 
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vides a setting where the interaction between different holes 
can be studied, an area that has yielded interesting results re- 
garding escape rates from circular'^ and diamond"" billiards 
with two holes, and very recently regarding transmission and 
reflection rates in stadium billiards^' and the bouncer modeP^. 
A general relation between escape and diffusion has been es- 
tablished by the escape rate theory of chaotic diffusion, which 
provides exact formulas expressing transport coefficients in 
terms of escape rates in spatially extended systems with ab- 
sorbing boundary conditions""'-'-"*'-''. 

Much research has gone into studying the parameter depen- 
dence of the diffusion coefficient in simple one-dimensional 
maps'^"'^. For low-dimensional, spatially periodic chaotic 
dynamical systems the diffusion coefficient is often found 
to be a fractal function of control parameters, exhibiting 
non-trivial fine scale structure even in apparently simple 
examples The source of this fractality is typically ex- 
plained in terms of topological instability under parameter 
variation of the underlying dynamics '"''^^. However there ex- 
ist systems that display these hallmarks of fractality but nev- 
ertheless have a linear diffusion coefficient"^. Therefore there 
is still work to be done explaining the phenomenon of fractal 
diffusion coefficients in one dimension, let alone attempting 
to answer questions about higher dimensional, more physi- 
cal systems like sawtooth maps^**, standard maps^'' or particle 
billiards^"'^' where analytical results are lacking, as are an- 
swers to basic questions about the structure of the diffusion 
coefficient. Previous work has focused on deriving and un- 
derstanding the diffusion coefficient under smooth variation 
of control parameters of the dynamics'^. In this setting the 
reduced, modulo one dynamics of a system will change with 
parameter variation''*. Here we switch focus and study a sys- 
tem where the reduced dynamics does not change^^'^^ -'^. The 
main aim of this paper is to see how the diffusion coefficient 
varies with the size and position of a hole in a paradigmatic 
'toy model'. In the process we find new relations between 
diffusion and escape rates, between diffusion and periodic or- 
bits, and we extend the small hole theory of Refs.2 and 3 to 
diffusion. What we learn here can hopefully be transferred to 
more physically realistic dynamical systems such as particle 
billiards "-'-\ 

In Sec. II we define the dynamical system that will be the 
main object of our study. It is a simple piecewise-linear 
chaotic map of the real line which is a deterministic realiza- 
tion of a random walk. It is constructed by copying and peri- 
odically lifting the 'Bernoulli shift' or 'doubling map mod- 
ulo one' over the whole real line^^'^^-^-. Modeling a coin 
tossing process in terms of deterministic chaos, the doubling 
map modulo one is one of the simplest dynamical systems 
displaying stochastic-like properties making it an indispens- 
able tool for studying statistical mechanics in the setting of 
dynamical systems In addition, we choose the doubling 
map so that we can compare with the results on escape rates 
from Refs.2 and 3 where the doubling map modulo one was 
also focused upon. Furthermore, the invariant measure of 
the doubling map modulo one is simply Lebesgue, which 
helps make it amenable to analysis with the method we will 
employ. The process of copying and periodically lifting a 



map is the classical way to study chaotic diffusion in one 
dimension'^ '''. However we do not introduce diffusion into 
the system through variation of a control parameter such as a 
shift or by varying the slope. Rather we 'dig holes' into the 
map that serve as intervals where points can be iterated to a 
neighboring interval in analogy with the work in Ref.2. 

We then analytically derive the diffusion coefficient as a 
function of the size and position of a hole in this system via 
the Taylor-Green-Kubo formula" '^, in terms of a functional 
recursion relation. There are various methods for analytically 
deriving diffusion coefficients'^'''''^*' but the method we use, 
developed in Refs.18, 25, 27, and 32, is the best suited to this 
setting. In Sec. Ill we look at the analytical formulas derived 
in Sec. II and find that the diffusion coefficient varies as the po- 
sition of a hole is varied, in analogy with results on the escape 
rate. We also find that the diffusion coefficient decreases non- 
monotonically as the size of a hole decreases, a result that is 
different to the escape rate. We explain this result in terms of 
the complicated forward and back scattering rules associated 
with even simple dynamical systems'-'''^'^'*'^^. 

By using the same approach, similar results are obtained 
for the diffusion coefficient of maps where the holes are not 
placed symmetrically, and where the map generating the mi- 
croscopic dynamics is non-symmetric. Following this, we 
consider the case of small hole size by deriving analytical 
expressions for the diffusion coefficient which capture the 
asymptotic regime. We find that the asymptotic regime is de- 
pendent upon the orbit structure of the limiting point in an 
escape region, a result which goes beyond a simple random 
walk approximation'^''^''^'^"*. We explain the results on po- 
sition dependence, non-monotonicity and asymptotic regimes 
by looking at the periodic orbit structure of the map. More- 
over, we build a periodic orbit expansion for small but finite 
holes giving a more intuitive insight of the above. In Sec. IV 
we numerically calculate the escape rate for the corresponding 
open system in order to compare with the structure of the dif- 
fusion coefficient. We summarize our findings and conclude 
our work in Sec.V. 



II. DERIVING THE DIFFUSION COEFFICIENT OF 
A DETERMINISTIC DYNAMICAL SYSTEM. 

The dynamical system that we will study is based on the 
doubling map modulo one, whose phase space is simply the 
unit interval. 

The tilde in Eq.(l) will be used throughout to signify a self- 
map. We turn Eq.(l) into a dynamical system that exhibits 
diffusion in two steps. Firstly, we dig two symmetric holes 
into i\7(x). Let < fli < 02 < 1/2 < 03 < 04 < 1, with 
04 = 1 — oi anda3 = 1— 02. For simplicity we let /i = a2— 04 
which is the size of a hole. We lift the map dynamics by 1 for 
a; G [ai, 02] and we lower the dynamics by 1 for x G [03, 04] 
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to create a map M{x) : [0, 1] [-1, 2], 



M{x) = < 



2x 
22; + 1 

2a; 
2x-l 
2x-2 
2x-l 



< X < ai 
ai < X < a2 
a2 < X < ^ 
^ < X < 

< X < ai 
a4 < X < 1 



(2) 



We label the intervals II ~ [01,02] and In ~ [03,04] for 
convenience. We call II and 1^ holes as they allow points to 
escape from the unit interval to a neighboring interval. 

Secondly, we periodically copy Af (x) over the entire real 
line with a lift of degree one such that. 



M{x + b) = M{x)+b, beZ, 



(3) 



so that M (x) : M ^ M. A uniform distribution of points on 
the unit interval will spread out when iterated under Eqs.(2,3). 
The diffusion coefficient D, is defined as the linear increase in 
the mean square displacement of a distribution of points and 
is given by the Einstein formula in one dimension as 



D = lim 



2n 



(4) 



where .t„ is the position of a point xq at time n which is given 
by M"{xo) in the system we consider. The angular brack- 
ets represent an average over a distribution of points. In the 
setting we consider, this distribution is the invariant density 
of the system p*{x) — 1, and the average we interpret as an 
integral. 



(...) = / ...p*ix)dx. 
Jo 



(5) 



Eq.(4) can be rewritten in terms of the velocity autocorrelation 
function of the system as the Taylor-Green-Kubo formula" '^, 

D = lim {voixHix))] - i {vo(xf) , (6) 

\k=0 / 

where Vk{x) = lxk+i\ — [xk\ gives the integer value of 
the displacement of a point xq at time k. Considering Eq.(2) 
Vk{x) takes the form. 



< Xk < ai 

1 Oi < Xk < 02 

Vk{x) = <( 02 < Xfc < 03 

-1 O3 < Xk < O4 

04 < Xfc < 1 



(7) 



The leading order term of Eq.(6), Drw, is simply equal to 
1 



Drw = 7; I ^'o{x) dx, 
^ Jo 

(04 - 03) + (02 - 01) 



= h. 



(8) 



Eq.(8) is the simple random walk result for diffusion that one 
obtains if higher order correlations are neglected'^ '^''^ -''*. In 



order to fully evaluate Eq.(6) we define a recursive function 

J"(a;) : [0, 1] ^ Z'«'25.2732^ 



r{x) = Y,Mx) 



k=0 



= Vo 

k=0 

= i>o(x) + J"-1(a7(x)). 



(9) 



We then define a cumulative function which integrates over 
Eq.(9) as in Eq.(6), 



T{x) = lim T\x) = / ,r{y)dy. 



(10) 



Due to the chaotic nature of the map M{x), J"(x) will be a 
very complicated step function for high values of n, hence in 
the limit n ^ 00 T{x) will be a fractal function exhibiting 
non-trivial fine scale structure"''^''^'^^. By combining Eq.(9) 
and Eq.(lO) we can solve T{x) as a functional recursion rela- 
tion. We use the conditions that r(0) = r(l) = and that 
the function T{x) is continuous to obtain 



hTi2x) 



Tix)=< ? 



T(2x) 
T{2x) 



x — ai 
02 - Ol 



^T{2x — 1) + 02 — Ol 
|T(2a; - 1) + 1 - a; - Ol 
jT(2a;- 1) 



< a: < Ol 

01 < a; < 02 

02 < a; < ^ 
i < a; < 03 

03 < a; < 04 

04 < .T < 1 



(11) 

Repeated application of the recurrence relation means we can 
solve Eq.(l 1) as an infinite sum, 

n 

T{x) = lim t{MHx)), (12) 



n— ^00 ^ — ^ 2 

fc=0 



where 



t{x) 




< X < Ol 

01 < X < 02 

02 < X < 03 

03 < 2^ < 04 
04 < X < 1 



(13) 



Eq.(6) can now be evaluated in terms of the functional recur- 
sion relation of Eq.(l 1) as 







= lim 




n— ^00 






( r- 


= lim 









(yj^ Mx)Y^^Vkix)dx^ ^0 



J"(x)dx- / J"(x)dx 



{x)dx 



r(o2) - r(ai) - r(o4) + T(a3) - h. 



(14) 



Finally, due to the condition that and Ir are symmetrically 
positioned, T(x) is a symmetric function. We can use this to 
simplify Eq.(14) to 



D = 2T(o2) - 2r(oi) - h. 



(15) 
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Eqs.(12,15) provide us with a very efficient way to evaluate 
the diffusion coefficient for any choice of position or size of 
II and Ir. For a more detailed discussion of this method see 
Refs.13, 18, 25, 27, and 32. We will evaluate Eq.(15) for a 
series of choices in the following section. 



III. ANALYZING THE DIFFUSION COEFFICIENT 

In this section we look at how the diffusion coefficient 
varies with the position of the holes and the asymptotic be- 
havior as the hole size goes to zero. 

A. Position dependence 

We evaluate Eq.(15) in a specialized setting where we re- 
strict II and to Markov intervals. That is, we choose the 
points fli, (i G {1,2,3,4}) to be dyadic rationals, i.e., for 
some fixed integer s > the points a.i are of the form 
with r £ Z and < r < 2". The points a.i will then be pre- 
images of 1/2 under the map M{x). As T(l/2) = 02 — ai, 
T{ai) can be evaluated with a finite sum rather than the infi- 
nite sum of Eq.(12). For each value of s there are 2""^ places 
to position an interval II of size 2^', with In being deter- 
mined by the symmetry condition. We can evaluate the diffu- 
sion coefficient at each of these choices via Eq.(15) and com- 
pare the results as the choices vary. For example, when s = 1 
there is only one choice for 1^, namely ai = 0, 02 = 1/2, ex- 
actly corresponding to a simple random walk, with 1/2 prob- 
ability of moving left or right at each step. The diffusion co- 
efficient for this system is well known to be 1/2, in agree- 
ment with the more general expressions given here, Eq.(8) for 
h = l/2andEq.(15), as 

D = 2T(l/2) - 2T(0) - (1/2) 

= 1/2. (16) 

For higher values of s the diffusion coefficient varies with the 
position of the holes, see Fig. 1. We see a step function that 
behaves increasingly erratically as the partition is refined and 
s is increased. We further note that the average of this step 
function can be calculated to be {Ds) = 2^* = h for a given 
s, which is the simple random walk solution of Eq.(8). 

The structure of the step functions in Fig. 1 can be explained 
in terms of the periodic orbits of the map M{x) which cor- 
respond to standing or running orbits of A/(a;)'-''''*'-'^. For 
example, if an image of II overlaps In, one will find a lot 
of backscattering in the system, i.e., points that escape the 
unit interval via II find themselves getting sent back via In 
(and vice versa). This has the result of decreasing the diffu- 
sion coefficient relative to the random walk solution derived 
in Eq.(8). In order to find intervals where this overlap occurs, 
we look for standing orbits by solving the simple equation 
MP{x) = X where M'^{x) = 1 — x for q < p. Due to the 
symmetry of the holes, the image of II containing the solu- 
tion of this equation will, after q iterations, overlap with In 
and backscattering will occur. The smaller values of p will 



correspond to values of x which give the most overlap and 
hence the most backscattering. For example, 

Af(a;) = l-a;, a; G [0,1/2], ^ a; = i. (17) 

Therefore if one places II so that x = 1/3 is in its interior, 
one will find the system has a relatively small diffusion coef- 
ficient due to the backscattering. This phenomenon caused by 
standing orbits'-''''^ is highlighted in Fig. 1. 

Alternatively, if the image of II overlaps with itself con- 
sistently then one will find a higher diffusion coefficient. 
This is due to the presence of running orbits or accelerator 
modes'^ ''* -'^ in such a system. In order to find such orbits, we 
solve the simple equation M'p{x) = x where M'^{x) 7^ 1 — x 
for q < p. For example p = 1 gives 

M{x)=x, a; g [0,1/2], =^ x = 0, (18) 

and we can see in Fig. 1 that when the hole contains the point 
X = one has a high diffusion coefficient relative to the sim- 
ple random walk result. When p = 2, 

M\x)=x, ^x = 0,i. (19) 

we can immediately throw the solution x = 1/3 away as this 
result corresponds to M(x) = 1 — x. However, forp = 3 

M\x)^x, ^x = 0,i (20) 

and again we see in Fig. 1 that these values correspond to rel- 
atively high diffusion coefficients when they are in the interior 
of II. This process of pinpointing standing and running peri- 
odic orbits can be continued for higher iterations with relative 
ease as we are dealing with a full shift map and there is no 
need to prune any solutions. This technique helps explain the 
increasingly complicated step function that one obtains as s is 
increased. 

At first sight, the step functions illustrated in Fig. 1 do not 
appear to contain much interesting structure. However, upon 
closer inspection we notice that every 'parent' hole of size 
2^* and associated diffusion coefficient Dg splits into two 
'child' holes of size 2"'^^+^^ and associated diffusion coef- 
ficients -Ds_(.i and D].^i respectively, such that 

D, = 2D'i+^ + 2Dl+^-2-\ (21) 

where superscripts 0, 1 correspond to left and right child hole 
respectively. To see this, one first needs to define cumulative 
functions T^{x) and T^(x) for the respective left and right 
child holes. Now since the cumulative functions are additive 
with respect to the holes we have that T{ai) ~ T'^{ai) + 
T^(ai) for i = 1 ... 4. Moreover, since the iterate of the par- 
ent hole endpoint (n always avoids both parent and child holes 
then T^{ai) = T^{ai). Finally, considering the midpoint 
ttm ~ {0.2 + 0.1)/ 2 of the parent hole which is also the right 
and left endpoint of the left and right child holes respectively, 
it follows from Eq.( 11) that TO (a„)-Ti(a„) = (a2-ai)/2. 
Eq.(21) follows after expanding in terms of r° and and 



4 



0.6 
0.5 
0.4 
0.3 
0.2 





(a) 









D 





0.18 
0.16 
0.14 
0.12 
0.1 
0.08 
0.06 H 
0.04 
0.02 








0.05 



0.04 



0.03- 



D 



0.02 



0.01 



0.25 

It 



0.125 0.25 0.375 

It 



0.3- 



0.2 



D 



0.1 



(c) 





0.5 

0.10- 









(b) 

















0.08 



0.06 



D 



0.04 



0.02 



0.5 



(e) 



0.125 



0.25 

It 



0.125 0.25 0.375 0.5 

It 



Id) 



0.375 



0.5 



0.125 0.25 0.375 

It 



0.5 




Figure 1. The diffusion coefficients: In these figures the diffusion coefficient D is illustrated for the doubling map M(a::) as a function of the 
position of the hole II of size 1/2". In (a), (b), (c), (d), (e) and (f) it is s = 2, 3, 4, 5, 6 and 12 respectively. D is given by the thick black lines 
whilst the holes are highlighted by the thin vertical lines. The thin horizontal lines are a guide to show the average value 1 /2* . The symbols in 
(f) refer to specific periodic orbits as discussed in the text. 
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substituting the above relations. Notice that recursive itera- 
tion of Eq.(21) n times gives an expression for Dg in terms of 
the 2" child diffusion coefficients 

I?, = (l-2")2-^+2" J2 ^'s+n. (22) 

J6{0,1}" 

where the sum runs over all 2" binary permutations of length 
n. Rearranging this we find 



E 

je{o,i}" 



- 2- 



9 — s — n 



(23) 



that is, the relative deviation of each diffusion coefficient from 
its mean is exactly additive. 

The above scaling and self-similarity structure (often con- 
sidered to be properties of fractal structures) of the step func- 
tions illustrated in Fig. 1 can be further investigated by defin- 
ing a set of continuous, cumulative functions which integrate 
over the step function, in the same way that T{x) integrates 
over the step function J"(a;). In order to define such a func- 
tion, $s (a;), for a given s, we firstly subtract the average diffu- 
sion coefficient {Ds) = 1/2'*, and integrate over the resulting 
step function. We then normalize this integral by multiplying 
it with 2''+^ so that it can be easily compared with other values 
of s. Let 



{D{y) 



dy, 



(24) 



where D{y) refers to the diffusion coefficient of the dyadic 
interval 1^ containing y. The solution to Eq.(24) is illustrated 
for several examples of s in Fig. 2. We see that as s increases, 
Eq.(24) becomes a fractal function exhibiting non-trivial fine 
scale structure and regions of scaling and self-similarity. This 
structure is symptomatic of the dense set of periodic orbits 
which exists in M{x). In the limit of s going to infinity, each 
periodic orbit makes the diffusion coefficient deviate from the 
average hence one obtains a dense step function. When this 
function is integrated over one sees a function that contains a 
dense set of maxima and minima, hence a fractal. 

Another interesting feature that we find in this system is 
that reducing the size of the holes can sometimes have no ef- 
fect on the diffusion coefficient. As we calculated in Eq.(16), 
if the holes are II — [0, 0.5] and Ir = [0.5, 1], the diffusion 
coefficient is equal to 0.5. However, we can reduce the hole so 
that II = [0, 0.25] and Ir = [0.75, 1] and the diffusion coef- 
ficient remains equal to 0.5 as illustrated in Fig. l.(a). We also 
see that if the holes are II = [0.25, 0.5] and Ir = [0.5, 0.75], 
the diffusion coefficient is equal to 0. This is due to a simple 
trapping mechanism in which no diffusion occurs. In addi- 
tion, reducing the size of the holes can result in an increase 
of the diffusion coefficient, i.e., the diffusion coefficient de- 
creases non-monotonically in some regions as the size of the 
holes is decreased. One can check this by comparing the fig- 
ures in Fig. I. This feature is due to the fact that increas- 
ing the size of the hole can in some cases introduce more 
backscattering into the system thereby reducing the diffusion 
coefficient: Consider the case where = [0.125,0.25] with 



diffusion coefficient equal to 1/16, Fig. l.(b), and compare 
with the case with smaller holes, II = [0.125,0.1875] but 
larger diffusion coefficient 5/64, Fig. l.(c). The dominant 
periodic orbit (the one with the lowest period) in the inter- 
val [0.125, 0.1875] is the orbit of the point 1/7 which corre- 
sponds to a running orbit as can be seen in Fig. l.(f). The 
dominant periodic orbit in [0. 1875, 0.25] meanwhile is the or- 
bit of the point 1/5 which corresponds to a standing orbit as 
can be seen in Fig. I.(f). The effect of decreasing the hole 
from [0.125,0.25] to [0.125,0.1875] is to remove the stand- 
ing orbit of 1/5 and the backscattering associated with it. 

One can also observe this phenomenon by looking at the 
fractal structure illustrated in Fig. 4.(b). While ai ~ 1/3 is 
fixed, the various maxima and minima that we see can be ex- 
plained by looking at the orbit of the point a2- We see that 
when 02 = 5/12 (corresponding to h = 1/12), the orbit of 
02 is a standing orbit and hence we see a striking minimum in 
the diffusion coefficient. If we reduce h so that h = 1/15 with 
02 = 2/5, the orbit of 02 is now a running orbit and we ob- 
serve a maximum in the diffusion coefficient. These points are 
highUghted in Fig. 4.(b). This explanation in terms of topolog- 
ical instability under parameter variation is discussed further 
inRefs.l3, 18,24, and 27. 



B. Non-symmetric models 

One can also consider the case where the holes are not 
placed symmetrically. One way of doing this is by setting 
ai = in Eq.(2) and letting Ir vary independently such that 

= ai < 02 < 1/2 < 03 < 04 < 1 and 04 — 03 = 
a2 = h = l/2•^ As before we assume Markov holes, i.e. 



03 = 1/2 + (i - l)h with i = 1,. 
function in this setting is given by 



iT"^(2x) 



+ X 

T"*(2x) + a2 
T""(2x-1)4 
T"^(2x - 1) - 
T"*(2x - 1) 



(2/1) . The cumulative 



< X < a2 



^2 



04 



a2 < X < ^ 
i < x < 03 
as < X < a4 
tti < X < I 



(25) 

where the superscript ns refers to the non-symmetrically 
placed holes. The diffusion coefficient is then given by 



D = T"'{a2) - T"'{ai) + Vias) - h , 



(26) 



cf. Eqs. (14,15). The consequence of using non-symmetrically 
placed holes can be seen in Figs. 3. (a), (b) and (c), in which 
the analytical Eq.(26) is evaluated for s = 4, 5 and 12 re- 
spectively and illustrated as a function of the position of Ir. 
We see that the step function structure is changed relative 
to Fig. l.(c), (d) and (f). In particular, in Fig. 3.(c) we ob- 
serve a very evenly distributed set of minima, these being the 
dyadic rationals which are in this case all the preimages of 
X = in and hence indicative of increased backscattering. 
In addition, the average value for the diffusion coefficient is 
(Dg) = 2(h — h?) rather than h due to the constant presence 
of the running orbit at x = 0. This follows from noticing 
that Ti(a4) — Tij^i{as) = —h, where Ti corresponds to a 
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X X 



Figure 2. (Color online) Cumulative integral function <&s(a;): In this figure the self similarity and scaling that one sees by integrating over 
the position dependent diffusion coefficient for the doubling map is illustrated. In (a) the structure is seen emerging as the hole size 1/2* is 
decreased. From (0.1, 0) upwards, s = 2 (grey), s = 5 (red), s = 8 (blue) and s = 20 (black). In (b) the region highlighted in (a) is blown up 
whilst the inset shows the highlighted region in (b) blown up in order to illustrate the self similarity and non-trivial fine-scale structure of the 
diffusion coefficient. 



cumulative function defined by 03 = 1/2 + (z — l)h with 
i = 1, . . . (2/1)^^. Here, the a4 for i is the same as 03 for i + 1 
and so their itinerary is the same except for the very first point 
where ti{a4) — ti^i{a3) = —h. So we have that the average 
over all I's is also equal to —h leaving (Dg) ~ T{a2) which 
can be calculated directly to be T{a2) = 2(/i — h^). 

Accordingly, one can now consider non-symmetric micro- 
scopic dynamics. One possibility for this is to replace the dou- 
bUng map modulo one in Eq.(l) by the tent map 



A(.t) : [0, 1] ^ [0, 1] , k{x) = 



2.T < a; < i 
i < .T < 1 



2x 



(27) 

Now the dynamics no longer commutes with the symmetry 
a; — > \ — X. We again dig two holes with < ai < a2 < 
l/2<a3<a4<l into Eq.(27) as before and periodically 
copy the resulting map over the real line with a lift of degree 
one. The cumulative function is given by 



\Ta{2x) 
Ta{2x)+x- 



ai 



Ta{2x) 



0.2 - «! 

-^Ta{2 - 2x) + a2 
-|Ta(2-22;) -x4 
-|Ta(2 - 2x) 



- ai 
04 



< a; < ai 
ai < X < a2 
a2 < X < ^ 

1 < X < 03 
< X < a4 

a4 < X < 1 



(28) 

where the subscript A indicates we are considering the cumu- 
lative function for the tent map. The diffusion coefficient is 
similarly given by 

D = TA{a2) - TA(ai) - rA(a4) + 1^(03) - h. (29) 

cf. Eq.(14). We can not use the symmetry condition to sim- 
plify Eq.(29) further as the cumulative function Eq.(28) is not 



a symmetric function. 

Interestingly, if the holes are placed symmetrically such 
that 04 = 1 — ai and 03 = 1 — 02, then D = h and is in- 
dependent of the hole position. This follows by noticing that 
TAias) = -l/2TA(2a2) + h andTAiai) = -l/2TA{2ai). 

If the holes are placed in a non-symmetric way such that 
^ ai < a2 < 1/2 < < < 1 and 04 — 03 = 02 = 
h = 1/2'*, then the periodic orbit dependent structure of the 
diffusion coefficient is re-established. This is illustrated as a 
function of Ij^ in Figs. 3.(d), (e) and (f) for s = 4, 5 and 12 
respectively. As before, {Ds) ~ Ta{0'2) = 2(/i — h'^). 

A recursive relation of the type of Eq.(21) can also be es- 
tablished for both doubling and tent maps with holes placed 
non-symmetrically, as described above. We find that 

2-s+i 



(30) 



where the difference in the last term's exponent is indicative of 
the constant running orbit at x = so that the average rescaled 
fluctuations of the child diffusion coefficients are intensified 
by a factor of two. 



C. Asymptotic behavior 

In this subsection we focus our discussion on the first case 
of symmetric holes in the symmetric map M{x) and comment 
on its generalizations to the non-symmetric cases towards the 
end. We will analyze the behavior of the diffusion coefficient 
for the Bernoulli shift as the hole size h = 02—01 goes to zero. 
By doing this, the hole will converge to a point which could 
be a running orbit, a standing orbit, or a non-periodic orbit. 
We derive equations which give the asymptotic behavior in all 
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Figure 3. Non- symmetric holes for a symmetric map and a non-symmetric map: In this figure the analytically derived diffusion coefficient 
D is illustrated for two different dynamical systems as a function of the position of the hole In, of size 1/2" (cf. Fig. 1) with II fixed at 0. 
We display the diffusion coefficient Eq.(26) of the doubling map Eq.(2) with s = 4, 5 and 12 in (a), (b) and (c) respectively and the diffusion 
coefficient Eq.(29) of the tent map Eq.(27) with s = 4, 5 and 12 in (d), (e) and (f) respectively. The horizontal lines give the average value 



three cases and use them to obtain the diffusion coefficient in 
terms of all periodic orbits in a hole of small but finite size. 
In order to do this, we first rewrite the (^(02) — T{ai)) term 
from Eq.(15) with Eq.(12) to 

T{a2)^T{ai) = 
Ji.- E ^ (^^' - i (fi' («i))) ■ (31) 

fc=0 

First consider the case that converges to a running orbit, 
that is, a periodic point of period p, which does not enter 
under forward iteration. In this case, from Eq.(13) we see that 
the only contributions to Eq.(31) come when k = lp,l ^ N 

Tia,)-na^)^limh(j2j^^ (h ^ 0) 

Evaluating Eq.(15) with Eq.(32) we get 

D{xp)^hj;^hf^±^\ (h^O), (33) 



where the superscript r denotes a running orbit. Now consider 
the case where II converges to a standing orbit, that is, a pe- 
riodic point Xp of period p, which enters Ir under forward 
iteration. Note that due to the symmetry of the holes this will 
always occur at time p/2 and hence standing orbits always 
have even periods. In this case we get a positive contribution 
to Eq.(31) when k = Ip, and a negative contribution when 
k = lp/2, 

T{a2) - T{a,) = lim h ( ^ ^ ) {h ^ 0) 

\/=0 ^ / 

^^{iTF^) ^^^^ 

In this case Eq.(15) evaluates as 

D{xp)^hj; = h[ —] (h^Q). (35) 

VI + 2 2 / 

The final case to consider is where II converges to a point 
which is non-periodic. In this setting the only contribution to 
Eq.(31) comes from the fc = term and therefore 

Dr^rh = h (h^O), (36) 
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which reproduces the simple random walk result. In summary, 
we have 



1-2-P 

2-P/2 



h p : 



: S 

n 



(37) 



Eq.(37) gives us a good explanation for the structure that we 
see in Fig. 1 with improved agreement for small holes (large 
s). As s is increased, the different asymptotic regimes can be 
seen in the step function. For example, when II is placed on 
a running orbit such as x = where p — 1, Eq.(33) tells us 
that D = 3h for small h. When II is placed on a standing 
orbit Uke x = 1/3 with p = 2, Eq.(35) tells us that D = h/3 
for small h. These deviations from the average value of h are 
observed in Fig. 1. 

A further consequence of Eq.(37) is the intriguing result 
that for small hole size, one can not rely on the simple ran- 
dom walk approximation for an accurate description of the 
diffusion coefficient'^''^'^^'^^. Rather, one must go beyond this 
theory and take into account the periodic orbit structure of the 
system, and in particular, the periodic orbits contained in the 
holes. The asymptotic regime that one obtains for small h will 
be dependent upon the type of point that the holes converge to. 
The authors are aware of only one other published result on a 
one-dimensional system in which the random walk approxi- 
mation theory is violated^^. In this case the phenomenon was 
explained in terms of ergodicity breaking, which is not the 
case here. 

We can now go beyond the small hole limit by combin- 
ing the above results with the parent-child hole relation of 
Eq.(21). For large n we have that 



Ds+n — 2 



(n -> oo) , 



(38) 



with p G {r, s, ri} depending on the limiting point of the hole 
as in Eq.(37). Hence, we may now express Dg in terms of all 
periodic orbits of period p which intersect the holes. 



{n 



(39) 



Note that a periodic orbit that intersects the parent hole more 
than once just gets added each time. Also, as discussed be- 
low, all periodic orbits are counted as running if they occur 
at the end of the interval. Eq.(39) suggests that the observed 
fluctuations of Dg from its average {Ds) = are due to the 
individual fluctuations of the infinitely many periodic orbits 
which intersect the holes. 

As expected from periodic orbit theory'"^, a very large num- 
ber of periodic orbits is needed to trace the hole accurately. 
However, if the periodic orbits are ordered appropriately the 
sum may be truncated to produce good approximations to 
Dg^^. Note that the optimal ordering (for fast convergence) 
of periodic orbits is by a modified version of the length of the 



orbit; J' 



pi 



J' 



^p, for 2pi « p2- In other words backscattering 
is much more dominant for orbits of equal period (see Fig. 4). 

We can further study the asymptotic behavior for the three 
different cases derived above and the finite hole size result by 



reducing II continuously. Fig. 4(a) illustrates these different 
regimes. 

When using Eq.(37), care needs to be taken when II con- 
verges to a point from the left or the right, i.e., it is not cen- 
tered on a point and reduced in size. In this case a boundary 
point of the hole, ai or a2, is kept fixed. If the boundary point 
is periodic, points near it in the interior of II miss Ip, and 
so it is always a running orbit. For example, as illustrated in 
Fig. 1 and Fig. 4. (a), x = 1/3 is a period two standing orbit 
when 1 /3 is in the interior of and the asymptotic regime for 
small h when J^, converges to 1/3 is given by Eq.(35) as h/3. 
However, if ai = 1/3 is fixed and h goes to zero, we must use 
Eq.(33) to evaluate the asymptotic regime as in this case 1/3 
is a running orbit. Eq.(33) tells us that the asymptotic regime 
is in fact bh/3. This additional topological subtlety that must 
be considered is illustrated in Fig. 4.(b). 

We close this subsection by considering small hole ap- 
proximations to the non-symmetric models considered in 
Sec. Ill B. In analogy to Eq.(37) for the doubling map with 
non-symmetric holes as in Eq.(25) and (26), we find through 
similar arguments that 



1 



1 



1 



(40) 



if //f is converging onto a running periodic orbit of period p 
(sinceT(a2) = 2hasa2 0). In contrast with the symmetric 
holes, here all backscattering is due to dyadic rationals. This 
is evident in Fig. 3. Hence, we find 



= 2 - 2^~ 



(41) 



if Ir is converging onto a dyadic rational of the form i/2'^' 
where i is some positive integer A non-periodic point gives 



j; = 2. 



(42) 



For the tent map with non-symmetric holes as considered 
in Sec. Ill B, Eqs.(40) and (42) carry over, while 



J" 



2 



(43) 



since dyadic rationals need one more forward iteration of the 
tent map to enter II ■ Finally, the periodic orbit sum of Eq.(39) 
modified by a factor of 2 also holds for the non-symmetric 
models considered here. 



IV. THE ESCAPE RATE 

It is interesting to note that an analytical relationship be- 
tween the escape rate of a spatially extended diffusive dy- 
namical system with absorbing regions and its diffusion co- 
efficient has been established by the escape rate theory of 
diffusion""'-' '*'^^. Motivated by Refs.2 where the compli- 
cated dependence of the escape rate on position and size of 
a hole has been studied, here we focus on the relationship be- 
tween the open map AI{x) on the unit interval with the sym- 
metric holes and serving as absorbing regions and the 
diffusion coefficient of the corresponding coupled, spatially 
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Figure 4. (Color online) The asymptotic regimes: In (a) the diffusion coefficient D{K) for holes centered on three different classes of points in 
M(x) is illustrated as a function of the hole size h (from top to bottom); a: = 1/3 a standing orbit (red), x = \/5/2 — 17/25 a non-periodic 
orbit (blue) and a; = 1/7 a running orbit (green) along with the different asymptotic regimes /i/3, h and 9/i/7 respectively shown by dashed 
lines. These asymptotic regimes correspond to the result of Eq.(37). In (b), the position of the left boundary of the hole a\ — 1/3 is fixed and 
h is again decreased continuously. We observe that a; = 1/3 becomes a running orbit when a critical point, and the asymptotic regime of 5/i/3 
(black dashed, top) illustrates this. The line h/3 (grey dashed, bottom) is what one would expect if x = 1/3 was contained in II yielding a 
standing orbit, and the random walk solution is given by the blue dashed line (middle). The two symbols (squares) identify parameter values 
where the right boundary point a2 of the hole generates a standing orbit, respectively a running orbit. 



extended system. That is, for calculating the escape rate any 
orbit that enters either of these intervals is removed from the 
system, and in this way points from an initial density escape, 
while for calculating the diffusion coefficient all points remain 
within the system by performing 'jumps' when hitting these 
intervals, as defined by the lift Eq.(3). An interesting question 
is to which extent the coupled diffusive 'jump dynamics' of 
the spatially extended system is already captured by the es- 
cape rate of the interval map that defines the unit cell of this 
lattice. 

The main result from Ref.2 concerning the escape rate is 
that escape will occur fastest through a hole whose minimal 
period is highest, or equivalently, the escape rate will be slow- 
est through the hole which has the smallest minimal period. 
By minimal period we mean the smallest period of all the pe- 
riodic points in a hole. 

In order to calculate the escape rate of our system we look 
at the transition matrix induced by the dynamics. The escape 
rate 7 can be evaluated via the largest eigenvalue ly of this 
transfer matrix'^''^''^'^"*-^^ 
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In v. 



(44) 



In Fig. 5 solutions to Eq.(44) are illustrated for s = 9 and 
compared with the diffusion coefficient in the corresponding 
extended system. Here we see that similar structures arise 
with deviations from the average occurring for both phenom- 
ena on the same intervals. In order to quantify these devia- 
tions we can compare Eq.(37) with Theorem 4.6. 1 from Ref.2, 



which generalized to Theorem 2.1 from Ref.3. It gives the es- 
cape rate for small hole size in the doubling map with one hole 
and can easily be generalized to escape through two holes as 
is the case here. This theorem states that the escape rate for 
small h with a running orbit (no iterate of the orbit reaches the 
second hole) is given by 



-^2 1- 



1 

2P 



(45) 



where Xp is the lowest period point in the hole with period p. 
For a standing orbit (the periodic orbit is in both holes) the 
period is effectively halved and we get 



2P/2 



{h 0) , 



(46) 



When the hole converges to a non-periodic point, the theorem 
states that the escape rate is given by 



h 



2 (/i ^ 0) 



(47) 



From Eq.(45) the relative deviation from the average escape 
rate (7) = 2h is given by 



lixp) - (7) 
for a running orbit and 

l{xp) - (7) = 



2h 
2P ' 

2h 
'2p7^' 



(48) 



(49) 
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Figure 5. (Color online) Comparing the diffusion coefficient with the 
escape rate: In this figure the diffusion coefficient for the doubling 
map is illustrated in black (bottom) alongside the escape rate of the 
corresponding open system in red (top) as a function of the hole II- 
The thin horizontal lines illustrate the average value to aid visual 
comparison of the fluctuations: (D) = 1/2'^ and (7) ~ 0.00393 
(3s. f.) ~ 1/2*. There is a clear relationship between the structure 
of these functions although intervals which give relatively high diffu- 
sion coefficients will give relatively low escape rates. The calculation 
has been performed for intervals II of size 1/2''. 



for a standing orbit. Similarly, the relative deviation from the 
average diffusion coefficient (D) = h, for a running orbit, can 
be obtained from Eq.(33) as 

D{xp) - (D) = (50) 

whilst for standing orbits, via Eq.(35), the relative deviation is 
given by 

2h 



D{xp) -{D) = - 



2P/2 + 1 ■ 



(51) 



Eqs.(48), (49), (50) and (51) help us explore the relationship 
between the diffusion coefficient of the extended system with 
the escape rate of the open system. An obvious difference 
is the absence of backscattering in the escape rate. However a 
more striking one is that the average escape rate does not equal 
the algebraic mean of all escape rates as for diffusion coeffi- 



cients. That is (7) 7^ 2=^^7=1 7s which is obvious from 



Fig. 5 but is also suggested by Eq.(45). The two symmetric 
holes are coupled differently for the escape problem and in a 
much more complicated way than as in Eq.(21) by involving 
the eigenvalues of 2""^ x 2''"^ transfer matrices. However, 
for small holes this coupling decays rapidly revealing the sim- 
ilarities which are seen in Fig. 5. We remark that while for 
non-symmetric holes as the ones considered in Sec. Ill B the 
diffusion coefficients may display qualitative differences (e.g. 
no position dependence), the corresponding escape rates differ 
only quantitatively. 



V. CONCLUSION 

The aim of this paper was to study the 'dependence of 
chaotic diffusion on the size and position of holes'. The 
answer was provided by analytically deriving the diffusion 
coefficient for both a symmetric and a non-symmetric one- 
dimensional piecewise linear map as a function of the size and 
position of a hole. We showed that for both maps the diffu- 
sion coefficient is a complicated function of the position and 
a non-monotonic function of the size of the holes, despite the 
fact that the underlying reduced dynamics is not changed, as 
is the case in previously studied models'^-^^-^^'^'. This finding 
implies that, surprisingly, making a hole smaller can increase 
the diffusion coefficient. These results we explained via the 
periodic orbit structure of the map and the ideas of running 
and standing periodic orbits. 

We furthermore found that the asymptotic regime that one 
obtains for small hole size is a function of the type of periodic 
orbit that the holes converge to. This is another important 
result, since it generalizes the standard uncorrected random 
walk approximation of simple stochastic processes. It implies 
that this random walk approximation may not always give ac- 
curate estimates for the diffusion coefficient of a chaotic dy- 
namical system. We have also obtained a new expansion for 
the diffusion coefficient of finite size holes in terms of periodic 
orbits and discussed their relative importance for the dynam- 
ics. In our setting, a periodic orbit can either be a running 
or a standing orbit. The presence of a standing orbit has the 
effect of reducing the diffusion coefficient relative to the aver- 
age value whilst the presence of a running orbit has the effect 
of increasing it relative to the average value. 

We finally numerically calculated the escape rate of the cor- 
responding open system and compared it with the diffusion 
coefficient thus relating diffusion and escape in a new man- 
ner We found that the diffusion coefficient and escape rate 
are both dependent upon the underlying periodic orbit struc- 
ture of the map, although differences arise which we explain 
as a difference in the coupling between holes. 

An interesting open question is whether there exists a 
parent-child scaling relation for the escape rate that is sim- 
ilar to the scaling relation for the diffusion coefficients 
Eqs.(21,30). Another interesting open question is to which ex- 
tent the above effects can be observed in computer simulations 
of diffusion in higher dimensional, more physically realistic 
systems such as suitably adapted periodic Lorentz gases'^''^ 
and related particle billiards'"''^". This should pave the way 
to design experiments where these effects might be observ- 
able, such as modified cold atom experiments on atom-optics 
bilUards^'^l 
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